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We shall revisit the conventional treatment of open quantum devices based on the Wigner- 
Function formalism. Our analysis will show that the artificial spatial separation between device 
active region and external reservoirs — properly defined within a semiclassical simulation scheme — 
is intrinsically incompatible with the non-local character of quantum mechanics. More specifically, 
by means of an exactly-solvable semiconductor model, we shall show that the application of the 
conventional boundary-condition scheme to the Wigner transport equation may produce highly 
non-physical results, like thermal injection of coherent state superpositions and boundary-driven 
negative probability distributions. 
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Current micro/nanoelectronics technology pushes de- 
vice dimensions toward limits where the traditional semi- 
classical Boltzmann theory can no longer be applied, 
and more rigorous quantum-transport approaches are re- 
quired 0]. However, in spite of the genuine quantum- 
mechanical nature of carrier dynamics in the active re- 
gion of typical nanostructured devices — like semicon- 
ductor superlattices and double-barrier structures — the 
overall behavior of such quantum systems is often the 
result of a non-trivial interplay between phase coher- 
ence and dissipation/dephasing [jj, the latter being also 
due to the presence of spatial boundaries It fol- 

lows that a proper treatment of such novel nanoscale 
devices requires a theoretical modelling able to properly 
account for both coherent and incoherent — i.e., phase- 
breaking — processes on the same footing. To this aim, 
a commonly used theoretical instrument is the so-called 
single-particle density matrix |J, |5j ; however, while the 
latter is the ideal tool for the description of ultrafast 
phenomena in infinitely-extended/periodically-repeated 
nanostructures Q , it cannot be directly applied to quan- 
tum systems with open boundaries, for which a real-space 
treatment is imperative. 

Such a real-space description is naturally provided by 
the Wigner-function formalism within this approach 
the statistical quantum state of the electronic subsystem 
is fully described in terms of the so-called Wigner func- 
tion, a function defined over the conventional phase-space 
as the Weyl- Wigner transform of the single-particle den- 
sity matrix 0, [nj. For the case of a one-dimensional 
problem, the equation describing the time evolution of 
the Wigner function — often referred to as the Wigner 
transport equation — is of the form: 



f(z,k) — / dz'dk' C(z, k; z' ,k')f(z' , k) , (1) 



d_ 

dt" 



where C{z, k; z', k') denotes the corresponding effective 



Liouville superoperator [lfj written in the Weyl- Wigner 
phase-space representation z, k. 

Different approaches for the study of quantum-tran- 
sport phenomena in semiconductor nanostructures based 
on the Wigner-function formalism have been proposed. 
On the one hand, starting from the pioneering work 
by Frensley a few groups ^3 have performed 

quantum-transport simulations based on a direct numer- 
ical solution of the Wigner equation in (JTJ via finite- 
difference approaches by imposing to f(z, k) the stan- 
dard boundary-condition scheme of the semiclassical de- 
vice modelling, also called U-scheme 0]. On the other 
hand, a generalization to systems with open bound- 
aries of the well-known Semiconductor Bloch Equations 
(SBE) 0, — the set of equations governing the time 
evolution of the single-particle density matrix — has been 
recently proposed |lCt uljj : such generalized SBE - 
obtained again via the U-scheme previously mentioned — 
describe the open nature of the problem via a bound- 
ary source term and a corresponding renormalization 
of the Liouville superoperator. In addition to the two 
alternative simulation strategies previously recalled - 
both based on effective treatments of relevant interac- 
tion mechanisms — , Jacoboni and co-workers have pro- 
posed a fully quantum-mechanical simulation scheme for 
the study of electron-phonon interaction based on the 
so-called "Wigner paths" E^. 

The generalized SBE approach previously recalled has 
been recently applied to prototypical semiconductor- 
based open systems. Preliminary results presented in 
|l3| seem to suggest an intrinsic limitation of the conven- 
tional Wigner function formalism in describing quantum- 
transport phenomena through systems with open bound- 
aries. On the other hand, no clear evidence of such 
limitations has been reported so far via Wigner-function 
simulations based on finite-difference treatments |3. H^. 
Aim of this Letter is to solve this apparent contradiction, 
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thus shedding light on the real limitations of the con- 
ventional Wigner-function picture applied to open-device 
modelling. Our analysis will show that the artificial spa- 
tial separation between device active region and external 
reservoirs is intrinsically incompatible with the non-local 
character of quantum mechanics. 

In order to gain more insight into this highly non- 
trivial problem, let us start considering the explicit form 
of the Wigner equation Q in stationary conditions and 
in the absence of any scattering mechanism: 

v(k) ^f(z, k)+ f dk'V(z, k - k')f(z, k') = 0, (2) 



where 



/-t-oc 
dz' 
-oo 



+00 e -2ikz' 



inh 



[V(z + z')-V(z-z')} (3) 



is the Weyl- Wigner superoperator corresponding to the 
device potential profile V(z), while v(k) denotes the elec- 
tron group velocity. Following the standard U-scheme, 
we shall now impose the desired spatial boundary con- 
ditions for / at the left (2 = — I) and right (z = +1) 
contacts, specifying the "incoming" electron distribution 
f b (k) = f(z b (k),k), where z b (k) = \ (6(-k) - 9(+k)) 
denote the left and right boundary coordinate corre- 
sponding, respectively, to positive and negative carrier 
wavevectors k (8 being the usual Heaviside step func- 
tion). By integrating Eq. J2J from the spatial boundary 
z b (k) to the current point z we get 



f(z,k) = f b (k) 



dz' f dk! — - — 1 —7—s -f(z',k') 

z b (k) J V{k) 



= fb(k)+ dz'dk'W{z,k;z',k')f{z',k') (4) 



with W(z,k;z',k') = -6(z< - z b {k))6{z - z') v{z ' ; k ( ~ k '\ 
In a compact notation we have: 



f = .fb + Wf, or / = 



1 



1 - W 



h 



(5) 



By expanding the above formal solution in powers of the 
interaction superoperator/propagator W — and therefore 
of the potential V — we get the well-known Neumann 
series: 



f = J2 wn h 



(6) 



n=0 



Let us now focus on the case of a symmetric po- 
tential profile (V(z) = V(—z)), which in turn corre- 
sponds to an antisymmetric potential superoperator, i.e., 
V(z,k) = —V(—z,k). Using this property together with 
the symmetric nature of our spatial boundaries, i.e., 
z b (k) = —z b (—k), it is possible to show that the interac- 
tion propagator W is also preserving the potential sym- 
metry 17]. This result is extremely important: it im- 
plies that for any symmetric potential profile and spatial 



boundaries the Neumann expansion in 10 gives always 
a Wigner function symmetric in the spatial coordinate: 
f(z,k) = f(—z,k). Therefore, in total agreement with 
the numerical results of the generalized SBE presented 
in |13| . the spatial charge density n{z) — J f(z,k)dk 
is always symmetric, no matter which is the shape of 
the injected carrier distribution f b (k). As anticipated, 
such symmetric behavior — which is an exact result of 
the treatment presented so far — has never been observed 
via finite-difference calculations. 

In order to solve this apparent contradiction, let us 
now focus on a particular choice of the device potential 
profile: V{z) = [1 + cos(2kz)]. The corresponding 
superoperator in © is simply given by: 

V(z, k) = ttV sm(2Kz)[5(k - k) - 5{k + k)} . (7) 

For this particular potential superoperator — characte- 
rized by a factorization/decoupling of position (z) and 
momentum (k) coordinates — it is possible to obtain the 
spatial charge distribution analytically: 



n(z) = / dkf b (k) iF 2 
with 



-; l + a(k), l-a(k); q(z) 



(8) 



q ( z ) = llUll sin [ K ( z _ !)] sin [ K ( z + !)] . (9) 



Here, m* is the electron effective mass, a{k) — k/ k, and 
1-F2 denotes the generalized hypergeometric function of 
type (1,2). 

In order to investigate the main features of the ana- 
lytical results obtained so far, let us start by consider- 
ing extremely simple spatial boundary conditions: a mo- 
noenergetic carrier injection from left only, i.e., f b (k) oc 
S(k — k ). Moreover, we choose a potential periodicity 
k = j corresponding to just one maximum within the 
device active region, so as to mimic the single-barrier 
device considered in 0. Figure ^ shows a compari- 
son between the analytical spatial charge distribution in 
© and the phenomenological result obtained via con- 
ventional scattering-state calculations. As we can see, 
the two curves differ significantly: while the phenomeno- 
logical charge distribution (dashed curve) — describing 
an extremely small tunnelling dynamics through the po- 
tential maximum — exhibits a clear and non-ambiguous 
asymmetric behavior, the analytical result of the Wigner- 
function approach (solid curve) is always symmetric, in 
total agreement with the analysis proposed in 01 ■ This 
implies that, also for an infinitely high potential barrier 
and for a monoenergetic injection from left only, carriers 
are "instantaneously" present also on the right part of 
the device, a typical non-local feature of our quantum- 
mechanical calculation. More important, within the 
Wigner-function treatment (solid curve) we deal with a 
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spatial charge/probability distribution with negative val- 
ues, which tells us that the analytical solution f(z, k) of 
the differential equation in 0) is not necessarily a Wigner 
function |Tsj . 

Let us finally come to potential discrepancies between 
exact and finite-difference results. Figure [21 shows a com- 
parison between the conventional finite-difference solu- 
tion of Eq. (0) proposed in 8] (dashed curve) and a 
numerical iterative solution of Eq. Q (solid curve); for 
both cases the same 80 x 80 phase-space discretization has 
been employed. As we can see, while the iterative solu- 
tion coincides with the analytical result (see solid curve 
in Fig. the finite-difference result comes out to be 
strongly asymmetric. A closer inspection (not reported 
here) reveals that such anomalous behavior is mainly as- 
cribed to the usual non-symmetric discretization (left or 
right derivative) of the kinetic/diffusion term in Eq. J2J, 
which may result in a fictitious decoherence/damping dy- 
namics |l9|. This seems to indicate that the more regu- 
lar — and physically-sound — results obtained so far via 
finite-difference calculations may be ascribed to such non- 
symmetric discretization procedure, which in turn tends 
to limit the highly non-physical features of the Wigner- 
function formalism applied to systems with open bound- 
aries. 

The analysis presented so far allows us to draw a few 
important conclusions. First of all, the results of the 
analytically solvable device model previously considered 
clearly show that the usual boundary-condition scheme 
— successfully applied to semiclassical device modelling — 
is intrinsically incompatible with the non-local nature of 
quantum mechanics. More specifically, from a strictly 
mathematical point of view it is true that the Wigner 
equation @ is a first-order equation in z parameterized 
by the wavevector k and, as such, it is sufficient to spec- 
ify the value of the Wigner function at the device bound- 
aries; on the other hand, from a physical point of view 
the separation between active region and external reser- 
voirs/contacts is only apparent, since the action of the 
potential superoperator V [see Eq. ©] is local in space, 
but its value inside the device depends on the properties 
of the potential v{z) both inside and outside the device 
active region [see Eq. We are forced to conclude 

that the application of the standard boundary-condition 
scheme to the Wigner equation in J5J is not physically 
justified, since it may provide solutions which are not 
Wigner functions, i.e., which do not correspond to the 
state of a quantum system. A clear and unambiguous 
proof of such non-physical outcomes are the negative val- 
ues of the electron probability distribution reported in 
Figs, n an d Generally speaking, what is intrinsically 
wrong in the usual Wigner-function treatment of open 
devices is the spatial separation between active region 
and reservoirs. This is similar to isolate a portion of a 
given energy spectrum, and try to treat such subset of 
energy levels as an independent subsystem. In contrast, 




position (nm) 

FIG. 1: Quantum-transport through a GaAs-based device 
active region (Z = 40 nm) with a cosine-like potential pro- 
file (V — 150 meV) sandwiched between its electrical con- 
tacts: comparison between the analytical charge distribution 
(solid curve) and the corresponding phenomenological result 
(dashed curve) for the case of a monoenergetic carrier injec- 
tion from left (E = 70meV). 
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FIG. 2: Finite-difference result (dashed curve) and a corre- 
sponding iterative numerical solution of Eq.[I](see text). 



in order to provide a rigorous description of open quan- 
tum devices it is imperative to introduce a clear sepa- 
ration between the degrees of freedom of the subsystem 
of interest (i.e., electrons primarily localized within the 
device active region) and those of the external reservoirs 
(i.e., carriers primarily localized within the electrical con- 
tacts). Given such separation, it is possible to obtain an 
effective transport equation for the subsystem of inter- 
est via a suitable statistical average over the coordinates 
of the "environment" . Following such reduction proce- 
dure, one gets again the Wigner equation where now 
the Liouville superoperator is the sum of a contribution 
describing the subsystem/device of interest plus an ef- 
fective superoperator describing the interaction of the 
device with external reservoirs: C — C dev + £, dev - res . 
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We stress that within the proposed alternative approach 
there is no need to impose any spatial boundary condi- 
tion: as for the case of a closed system the spatial domain 
is again — oo < z < +00. The only boundary condi- 
tion is the Wigner function j[z, k) at the initial time t ; 
it follows that — contrary to the case of spatial bound- 
ary conditions — the homogeneous Liouville dynamics in 
JU provides/maintains a "good" Wigner function at any 
later time t, thus preventing from non-physical behaviors, 
like those reported in Fig. 

We finally stress that the alternative simulation scheme 
described so far may be concretely realized following the 
prescription recently proposed in j^: the basic idea is 
to replace the usual open description of quantum devices 
based on spatial boundaries (i.e., source/loss terms) with 
a closed-system treatment, where the interaction of the 
device electrons with the external reservoirs is simply de- 
scribed in terms of effective scattering rates acting on the 
device electrons only. 
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